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1.  Productivity  Measures 


•  Number  of  book  published:  1 

•  Number  of  refereed  journal  articles:  14 

•  Number  of  articles  in  Conference  Proceedings:  12 

•  Number  of  refereed  papers  submitted/to  be  submitted:  5 

•  Number  of  graduate  students  supported:  4 

2.  Summary  of  Research  Tasks 

•  Absorptivity  of  Laser  Radiation  by  Metal  Powders 

-  Temperature  and  wavelength-dependent  spectral  absorptivities  for  some 
important  metallic  materials  in  the  infrared  were  studied. 

•  Modeling  of  SLS  of  Two-Component  Metal  Powders 

-  Analytical  models  of  melting  and  resolidification  of  powder  bed  under  constant  or 
temporal  Gaussian  heat  source  were  developed. 

—  Two-dimensional  models  for  melting  accompanied  by  shrinkage  and 
resolidification  were  developed  and  parametric  studies  were  performed. 

—  Three-dimensional  models  for  SLS  of  two-component  metal  powders  were 
developed  by  accounting  for  melting,  shrinkage,  capillary  flow,  and 
resolidification. 

•  Modeling  of  SLS  of  Single-Component  Metal  Powders 

-  Analytical  model  of  melting  of  powder  bed  under  constant  heat  source  were 
developed. 

-  Three-dimensional  model  for  SLS  of  single-component  metal  powders  were 
developed  by  considering  partial  melting,  shrinkage,  capillary  flow,  and 
resolidification. 

-  Three-dimensional  model  for  Selective  Laser  Powder  Remelting  (SLPR)  of 
single-component  metal  powders  were  developed  by  solving  melting  with 
shrinkage,  resolidification,  as  well  as  convection  in  the  liquid  pool. 

•  Modeling  of  Post-Processing  of  the  Laser  Sintered  Porous  Part 

-  A  model  for  liquid  metal  in  filtration  and  solidification  during  post-processing  of 
laser  sintered  porous  part  is  developed. 

3.  Detailed  Summary  of  Technical  Progress 

The  project  involves  fundamental  modeling  of  the  laser  beam-material  interactions  associated 
with  Selective  Laser  Sintering  (SLS)  of  single-  and  two-component  metal  powders,  as  well  as 
post-processing  of  the  sintered  parts  with  liquid  metal  infiltration. 

3.1  Absorptivity  of  Laser  Radiation  by  Metal  Powders 

Any  model  for  laser  processing  of  materials  must  have  a  complete  description  of  the  coupling 
between  the  laser  source  and  the  material.  The  coupling  is  defined  by  the  spectral  absorptivity  of 
the  material  at  the  wavelength  of  operation;  hence  the  normal  spectral  absorptivity  is  a  critical 
parameter  of  interest.  Optical  properties  of  bulk  metals  are  typically  functions  of  wavelength, 
temperature,  surface  geometry  (roughness),  incident  intensity,  and  physical  atomic  structural  and 
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electrical  properties  of  the  material.  Since  absorptivity  is  dependent  on  temperature,  it  is  also 
important  to  consider  the  absorptivity  change  in  the  modeling  of  laser  heating. 

Temperature  and  wavelength-dependent  spectral  absorptivities  for  some  important  metallic 
materials  in  the  infrared  were  studied  by  Boyden  and  Zhang  [J3,  C7]  .  For  the  application  of 
laser  material  processing,  temperature-dependent  values  are  calculated  based  on  available 
experimental  data.  Figure  1  shows  the  comparison  of  the  predicted  absorptivity  of  AISI  304 
stainless  steel  as  a  function  of  wavelength  with  available  experimental  data  in  the  literature.  The 
calculated  absorptivity  of  the  major  component 
in  AISI  304  -  Fe  is  also  shown  in  Fig.  1.  The 
absorptivity  of  AISI  304  is  much  higher  than 
that  of  pure  iron  metal.  It  can  be  seen  that  our 
calculated  results  agreed  very  well  with  the 
experimental  data.  The  optical  constants  and 
absorptivity  of  other  element  (Al,  Cu,  and  Ni) 
and  alloys  (Inconel,  A-110-AT,  Al  2024,  Al 
7075,  and  TiB120  VAC)  are  also  calculated 
based  on  the  Drude  type  model.  The  predicted 
absorptivities  of  the  pure  metals  at  10.6  pm 
agreed  very  well  with  the  experimental  results 
except  for  the  transition  metals  due  to  interband 
absorption.  Agreement  of  alloy  and  element 
absorptivity  calculated  values  and  the 
experimental  data  is  good  at  10.6  pm  but  not  at 
1 .06  pm.  The  higher  absorption  at  1 .06  pm  can 
be  contributed  to  parallel  band  absorption, 
which  is  not  included  in  the  Drude  absorption 
theory. 

3.2  Modeling  of  SLS  of  Two-Component  Metal  Powders 

Analytical  models  of  melting  and  resolidification  of  powder  bed. 

Melting  and  resolidification  in  SLS  of  two-component  metal  powders  in  SLS  processes  were 
systematically  investigated  by  the  PI  and  students.  Melting  of  a  subcooled  powder  bed  with  the 
finite  thickness  that  contains  a  mixture  of  two  metal  powders  with  significantly  different  melting 
points  is  investigated  analytically  [J6,  Cl].  Shrinkage  induced  by  melting  is  taken  into  account  in 
the  physical  model.  The  effects  of  porosity,  Stefan  number,  and  subcooling  on  the  surface 
temperature  and  solid-liquid  interface  are  also  investigated. 

In  order  to  discover  the  advantages  of  utilizing  a  pulsed  laser  in  a  SLS  process,  melting  and 
resolidification  of  a  subcooled  two-component  metal  powder  bed  with  a  temporal  Gaussian  heat 
flux  was  investigated  analytically  by  Konrad  et  al.  [J1 1].  The  laser-powder  bed  interaction  can 
be  divided  into  three  stages:  (1)  preheating,  (2)  melting  with  shrinkage,  and  (3)  resolidification; 
their  solutions  were  obtained  by  using  an  integral  approximate  solution.  The  results  show  that  an 
increase  in  heat  source  intensity  or  powder  bed  porosity  will  result  in  an  increase  of  the  melt  pool 
depth,  surface  temperature,  and  overall  processing  time. 


Wavelength  (jam) 


Figure  1  Wavelength-dependent  absorptivity 
for  AISI  304  stainless  steel  [J3] 


’  See  lists  of  journal  and  conference  papers. 
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Two-Dimensional  Modeling  of  SLS  of  Two-Component  Metal  Powders 

SLS  of  the  first  layer  is  modeled  as  melting  and  resolidification  of  a  metal  powder  layer 
subject  to  a  moving  heat  source  on  top  while  the  bottom  is  adiabatic  [J 1 ,  C3].  SLS  of  the 
consecutive  layer  is  modeled  as  melting  and  resolidification  of  a  metal  powder  layer  on  top  of 
the  existing  multiple  resolidified  layers  [Jl,  C2].  The  results  indicate  that  the  thicknesses  of  the 
loose  metal  powder  layer,  the  moving  heat  source  intensity  and  the  scanning  velocity  have 
significant  effects  on  the  sintering  process  in  both  the  first  layer  and  each  subsequent  layer.  A 
parametric  study  is  performed  and  the  best  combination  of  processing  parameters  is 
recommended. 

Three-Dimensional  Modeling  of  SLS  of  Two-Component  Metal  Powders 

The  highest  level  of  our  efforts  is  represented  by  modeling  of  three-dimensional  SLS  process 
with  both  shrinkage  and  liquid  flow  considered  [J4-5,  J9-10,  C6,  C9-11].  Laser  sintering  of  a 
two-component  metal  powder  layer  under  a  moving  Gaussian  laser  beam  was  investigated 
numerically  by  Chen  and  Zhang  [J4,  C6,  CIO],  The  laser  induced  melting  with  shrinkage,  and 
resolidification  of  the  metal  powder  layer  are  modeled  using  a  temperature  transforming  model. 
The  liquid  flow  of  the  melted  low  melting  point  metal  driven  by  capillary  and  gravitational 
forces  is  also  taken  into  account.  Simulations  were  conducted  using  a  mixture  of  40%  nickel 
braze  powder  and  60%  AISI  1018  carbon  steel  powder  by  volume.  To  validate  the  simulation 
results,  the  results  of  numerical  solution  are  compared  with  the  experimental  results  as  shown  in 
Fig.  2.  The  black  area  in  the  micrograph  is  a  local  void  and  the  light  area  is  sintered  metal.  It  can 
be  seen  that  the  actual  and  predicted  shapes  of  the  Heat  Affected  Zone  (HAZ)  are  similar,  but 
large  degree  of  local  porosity  exists  in  the  experimentally  observed  HAZ;  this  indicates  that  the 
HAZ  in  the  experiments  are  not  fully  densified  as  we  assumed.  Therefore,  a  partial  shrinkage 
model  was  developed  to  allow  part  of  the  gas  remains  in  the  HAZ  [J5].  The  results  showed  that 
the  predicted  results  agreed  very  well  with  experiments  if  it  is  assumed  that  the  volume  fraction 
of  the  gas  in  the  HAZ  is  0.2. 
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Figure  2  Comparison  of  numerical  and  experimental  solutions  [J4] 

In  reality,  SLS  is  a  layer-by-layer  process  by  which  the  sintering  process  occurs  in  a  fresh 
loose  powder  layer  on  top  of  multiple  sintered  layers.  Numerical  solution  of  three-dimensional 
laser  sintering  of  a  two-component  metal  powder  layer  on  the  top  of  multiple  sintered  layers  was 
investigated  for  the  case  of  single-line  scanning  [J9,  C9,  CIO]  and  multiple-line  scanning  [J10, 
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C 1 0- 11].  The  modeling  of  single-line  scanning  addresses  laser  scanning  in  the  middle  of  the 
loose  powder  on  top  of  the  existing  sintered  layers,  and  the  multiple  line  scanning  is  modeled  as 
sintering  at  the  boundary  of  loose  powder  layer  and  sintered  layer  as  shown  in  Fig.  3;  these 
treatment  are  reasonable  because  the  diameter  of  the  laser  beam  is  much  smaller  than  the  sizes  of 
the  powder  bed  in  the  horizontal  direction,  and  the  effect  of  laser  heating  is  limited  in  the  region 
very  close  to  the  laser  spot. 

Figure  4  shows  the  simulated  results  for  the  single  line  scanning  -  which  results  in 
axisymmetric  HAZ  about  Y=0  and  only  half  (Y>0)  of  it  is  shown  -  and  multiple-line  laser 
scanning  -  the  HAZ  is  no  longer  axisymmetric.  The  volume  fractions  of  the  gas  remaining  in  the 
HAZ  for  both  cases  were  <pg  t  =  0.2 ,  which  is  reduced  from  0.42  in  the  loose  powder. 


O  Low  Melting  Point  Powder 
•  High  Melting  Point  Powder 


Figure  3  Physical  model  of  SLS  by  multiple-line  laser  scanning  [J 10] 


(a)  Single-line  laser  scanning  (axisymmetric  (b)  Multiple-line  laser  scanning  (asymmetric 
about  Y=0)  [J9]  about  Y=0)  [J  1 0] 

Figure  4  Three-dimensional  shape  of  the  HAZ  ( (pg  ,  =  0.2  ) 
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3.3  Modeling  of  SLS  of  Single-Component  Metal  Powders 

Analytical  models  of  partial  melting  of  powder  bed 

The  final  product  of  SLS  of  the  two-component  powder  exhibits  the  mechanical  properties 
and  characteristics  of  their  weaker  component;  therefore,  the  interest  in  the  production  of 
metallic  objects  using  single  component  powders  has  been  increased  in  the  recent  years.  One 
approach  is  to  employ  partial  melting,  in  which  only  the  surfaces  of  particles  are  molten  and  the 
powders  are  sintered  by  binding  the  non-melted  solid  cores.  Since  SLS  is  a  very  rapid  process, 
the  liquid  layer  and  solid  core  of  a  partially  molten  powder  particle  are  not  at  thermal  equilibrium 
and  have  different  temperatures.  The  mean  temperature  of  a  partially  molten  grain  can  be  either 
below  or  above  the  melting  point  depending  on  the  degree  of  melting.  Thus,  the  melting  process 
can  be  assumed  to  occur  in  a  temperature  range  and  the  degree  of  melting  is  function  of 
temperature  in  the  range.  Partial  melting  of  a  single-component  metal  powder  with  constant  heat 
flux  was  investigated  analytically  by  Xiao  and  Zhang  [J13,  C5].  The  partial  melting  model  was 
also  extended  to  melting  of  alloy  by  including  a  completely-melted  liquid  layer  on  top  of  the 
mushy  zone  [J 1 4], 

Three-dimensional  model  for  SLS  via  partial  melting 

Xiao  and  Zhang  [J7,  C8]  numerically  solved  three-dimensional  partial  melting  of  a  single¬ 
component  metal  powder  bed  subject  to  a  moving  Gaussian  laser  beam.  The  liquid  velocities 
induced  by  capillarity  and  gravitational  forces  and  solid  velocities  caused  by  shrinkage  during 
the  melting  process  were  taken  into  account.  Numerical  calculations  were  performed  for 
sintering  of  AIS1  304  stainless  steel  powder.  In  order  to  minimize  the  balling  effect  and 
geometric  distortion  (curling  or  delaminating),  a  pre-heated  powder  bed  with  initial  temperature 
of  380°C  was  used  in  the  simulation.  Figure  5  illustrates  the  surface  temperature  distribution  of 
the  powder  bed.  It  can  be  seen  that  the  peak  temperature  at  the  powder  bed  surface  is  near  the 
trailing  edge  of  the  laser  beam  rather  than  at  the  center  of  the  laser  beam  due  to  motion  of  the 
laser  beam.  Figure  6  shows  the  three-dimensional  shape  of  the  powder  bed  surface,  molten  pool 
and  HAZ  at  the  same  condition  of  Fig.  5.  The  shrinkage  at  the  track  of  the  center  of  the  laser 
beam  is  most  significant  and  a  groove  is  formed  on  the  surface  of  the  powder  bed. 


T 


Figure  5  Temperature  distribution  on  the 
powder  bed  surface  (Nt  =  0.\9,Ub  =  0.12)  [J7] 


Figure  6  Three-dimensional  shape  of  the 
HAZ  (N,  =0.19, Uh  =  0.12)  [J7] 
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Three-dimensional  models  for  Selective  Laser  Powder  Remelting  (SLPR) 

Modeling  of  SLS  of  single-component  metal  powder  trough  complete  melting  -  an  approach 
to  fabricate  fully-densified  single  component  part  -  is  also  performed.  Liquid  flow  occurs  in  a 
non-porous  liquid  pool,  and  the  effect  of  liquid  flow  driven  by  the  surface  tension  gradient  and 
gravitational  force  on  the  SLS  process  is  much  stronger.  The  Darcy’s  law  that  was  used  to  obtain 
the  liquid  velocity  for  the  case  of  partial  melting  is  no  longer  applicable  and  the  convection  in  the 
liquid  pool  must  be  described  by  the  complete  three-dimensional  Navier-Stokes  equations.  A 
three-dimensional  model  describing  melting  and  resolidification  of  direct  metal  laser  sintering 
process  under  the  irradiation  of  a  moving  Gaussian  laser  beam  is  developed  by  Xiao  and  Zhang 
[J12].  We  also  developed  a  three  dimensional  model  describing  melting  and  resolidification  of 
direct  metal  laser  sintering  of  loose  powders  on  top  of  sintered  layers  with  both  single-  [Cl  2,  SI] 
and  multiple-line  laser  scanning  [S3].  Effects  of  shrinkage  and  natural  convection  driven  by  the 
surface  tension  and  buoyancy  force  are  taken  into  account.  The  energy  equation  is  formulated 
using  temperature  transforming  model  and  solved  by  the  finite  volume  method.  Temperature 
distribution  and  velocity  field  are  investigated.  The  results  show  that  increasing  initial  porosity  of 
the  powder  bed  enlarges  the  depth  of  the  melt/solid  interface  and  laser  intensity  has  great 
influence  in  both  depth  and  width  of  the  liquid  pool. 


Figure  7  Three-dimensional  shape  of  the  HAZ  ( As  =  0.4 ,  Ub  =  0.02 ,  £  =  0.5  ,  N  =  5 )  [C 1 2] 


(a)  top  view  (b)  longitudinal  view  at  y  =0  (c)  cross-sectional  view  at  x  =0 


Figure  8  Dimensionless  velocity  vector  plot  ( As  =  0.4 ,  Uh  =  0.02 ,  s  =  0.5 ,  N  =  5 )  [C 1 2] 


8 


Figure  7  shows  the  three-dimensional  shapes  of  the  powder  bed  surface,  melt  pool  and  HAZ 
for  single-line  scanning  with  five  sintered  layers  underneath  (N  =  5)  and  an  initial  porosity  of 
0.5  (f  =  0.5).  The  depths  of  surface  and  HAZ  decrease  with  the  increasing  X  in  the  laser 
scanning  direction  Y.  Figure  8(a)-(c)  show  the  velocity  vectors  in  the  liquid  pool  plotted  in  three 
different  views.  Since  the  surface  tension  is  a  decreasing  function  of  temperature,  i.e., 
dy/dT  <  0 ,  the  higher  surface  tension  of  the  cooler  liquid  metal  near  the  edge  of  the  liquid  pool 
tends  to  pull  the  liquid  metal  away  from  the  center  of  the  liquid  pool,  where  the  liquid  metal  is 
hotter  and  the  surface  tension  is  lower.  Therefore,  fluid  flow  on  the  surface  of  liquid  pool  is 
radially  outward  as  can  be  seen  in  Fig.  8(a).  The  liquid  metal  flow  is  also  driven  by  buoyancy 
force  as  illustrated  in  Fig.  8(b)  and  8(c).  The  hotter  liquid  metal  near  the  central  region  of  the 
molten  pool  flows  up  to  the  surface,  while  the  cooler  liquid  metal  near  the  pool  boundary  sinks 
along  the  melt/solid  interface  to  the  bottom  of  the  pool. 

Figure  9(a)-(c)  present  the  temperature  contour  in  the  powder  bed  plotted  in  three  different 
views  (N  =  5).  The  Marangoni  convection  (radially  outward)  at  the  top  surface  and  the 
shrinkage  phenomenon  resulted  in  a  significant  amount  of  heat  flows  from  the  hotter  region  to 
the  colder  region  especially  in  the  z  direction,  which  in  turn  result  in  a  wider  and  deeper  melt 
pool.  Another  observation  is  that  the  isotherms  near  the  melting  front  are  more  closely  spaced 
compared  with  those  far  away  from  the  melt/solid  interface. 


(a)  top  view  (b)  longitudinal  view  at  y  =0  (c)  cross-sectional  view  at  x=0 

Figure  9  Dimensionless  temperature  contour  ( =  0.4 ,  Ub  =  0.02 ,  e  =  0.5 ,  N  =  5 )  [C 1 2] 

3.4  Modeling  of  Post-Processing  of  the  Laser  Sintered  Porous  Part 

The  parts  produced  by  SLS  with  single  or  multiple  components  powders  are  usually  not  fully 
densified  and  have  porous  structure.  In  order  to  produce  fully  densified  part,  post-processing  is 
necessary.  Compared  to  sintering  and  HIP  processes,  the  advantage  of  infiltration  is  that  the  full 
density  can  be  achieved  without  shrinkage  in  the  post-processing.  The  additional  advantages  of 
post-processing  with  infiltration  include:  it  is  relative  inexpensive,  and  tooling  is  similar  to 
casting  process.  Infiltration  is  a  process  where  a  liquid  metal  is  drawn  into  the  pores  of  a  solid 
(part  produced  by  SLS  of  metal  powder)  by  capillary  forces.  The  liquid,  as  it  advances  through 
the  solid,  displaces  vapor  from  the  pores  and  leaves  behind  a  relatively  dense  structure.  The  rate 
of  infiltration  is  related  to  the  viscosity  and  surface  tension  of  the  liquid,  and  the  pore  size  of  the 
SLS  parts. 

In  order  to  understand  the  liquid  flow  associated  with  post-processing,  a  numerical  study  of 
transient  fluid  flow  and  heat  transfer  in  a  porous  medium  with  partial  heating  and  evaporation  on 
the  upper  surface  is  performed  [J8].  The  dependence  of  saturation  temperature  on  the  pressure 
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was  accounted  for  by  using  Clausius-Clapeyron  equation.  The  results  showed  that  as  time  passes 
the  magnitude  of  velocity  increases  until  the  process  reached  steady  state.  Solidification  of  liquid 
copper  infiltrated  in  steel  skeleton  was  modeled  using  temperature  transforming  model  in 
conjunction  with  Brinkmann  model  for  liquid  metal  flow  in  porous  media.  The  model  was 
validated  with  the  experimental  data  available  in  the  literatures  and  parametric  study  was 
performed  [S4].  Infiltration  and  solidification  of  liquid  metal  into  laser  sintered  porous  structure 
is  modeled  using  Volume  of  Fluid  (VOF)  and  temperature  transforming  model  [S5].  The  Ph.D. 
student  (Piyasak  Damronglerd)  who  is  working  on  the  modeling  of  post-processing  is  scheduled 
to  graduate  in  Spring  2007  and  is  in  the  final  stage  of  completion  of  his  dissertation. 
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